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Abstract. We present a detailed study of the viability of general vector-tensor 
theories of gravity in the presence of an arbitrary temporal background vector field. 
We find that there are six different classes of theories which are indistinguishable from 
General Relativity by means of local gravity experiments. We study the propagation 
speeds of scalar, vector and tensor perturbations and obtain the conditions for classical 
stability of those models. We compute the energy density of the different modes and 
find the conditions for the absence of ghosts in the quantum theory. We conclude that 
the only theories which can pass all the viability conditions for arbitrary values of the 
background vector field are not only those of the pure Maxwell type, but also Maxwell 
theories supplemented with a (Lorentz type) gauge fixing term. 
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1. Introduction 

Although the interest in alternative gravity theories has been present over the years [HE], 
it has been recently renewed by the unexpected observations made in cosmological and 
astrophysical contexts. Indeed, the difficulties found by the Standard Model (SM) of 
elementary particles and General Relativity (GR) in order to explain the nature of dark 
matter and dark energy suggest either the need for new physics beyond the SM or the 
modification of GR at very large scales. 

Among such modifications, we find the so called vector-tensor theories of gravity, 
proposed in the early 70's [Sj H], and where, in addition to the metric tensor, a 
dynamical vector field is introduced in such a way that it only couples to matter through 
gravitational interactions. 

After some initial proposals, these theories were practically forgotten due to 
consistency problems and only particular classes of models with fixed norm vector or 
potential terms have been considered in the context of violations of Lorentz invariance 
(see |E] and references therein). However, very recently, it has been shown [HIT] that 
unconstrained vector-tensor theories without potential terms are excellent candidates 
for dark energy. 

The main three difficulties found in this type of theories were the following [H [9]: 

a) Inconsistencies with local gravity tests: in addition to modifications in the 
static parametrized post-Newtonian (PPN) parameters (7, /?), these theories typically 
predict preferred frame effects, parametrized by (01,02). Such effects are produced 
by the motion with respect to the privileged frame selected by the presence of the 
background vector field. The increasing observational precision in Earth and Solar 
System experiments seems to exclude this type of effects [2]. 

b) Classical instabilities: they are generated by the presence of modes with negative 
propagation speed squared. This gradient instabilities give rise to exponentially growing 
perturbations which make the theory unphysical. 

c) Quantum instabilities (ghosts): they typically appear when the energy of some 
perturbation modes becomes negative. In such a case it has been shown that the vacuum 
state of the corresponding quantum theory is unstable [10]. 

It has been also argued [TTj that the presence of faster than light propagation modes 
could give rise to causality inconsistencies although this conclusion is far from clear (see 
for instance [I2]) and even some authors claim that consistency requires the presence 
of such modes [9]. For that reason in this work we will not include it as a viability 
requirement. 

Although it is usually claimed that vector-tensor theories are plagued by the 
mentioned consistency problems, the actual situation is that none of them has been 
studied in detail, and they are only based on qualitative arguments which, in many 
cases, turn out to be inappropriate. In this work we present a detailed analysis of the 
viability of unconstrained vector-tensor theories and obtain the precise conditions in 
order to satisfy the above requirements. 
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2. Local gravity constraints 

Without any other restriction, apart from having second order hnear equations of 
motion, the most general action for a vector-tensor theory without potential terms 
can be written as follows [2]: 



with uj, a, T, 6 dimensionless parameters and F^,^ = df^Ai, — di^A^. Notice that the 
term including the Ricci tensor can be rewritten as V ^A^^V ^A'^ — ^ fj.A'^V ^A^^ . These 
two terms do not appear independently in the most general action because they can be 
recasted in the terms written in ([1]). 

Since GR agrees with Solar System experiments with high precision, it would be 
desirable that a viable vector-tensor theory had the same set of PPN parameters as GR, 
i.e. 7 = /? = 1 and cti = 02 = 0. Notice that the rest of PPN parameters vanishes 
identically for models described by i^. The PPN parameters for general vector-tensor 
theories have been calculated in [2] assuming the existence of a constant background 
time-like vector field: 



and can be found in the Appendix. When we impose that such parameters agree with 
those of GR for any value of A, we obtain two sets of compatible models, according to 
their behavior in flat space-time, both with u = 0: 

• Gauge non-invariant models. These models have r 7^ and the corresponding term 
V fj^Ay'V^A'^ breaks the U{1) gauge invariance in Minkowski space-time. The three 
possibilities we obtain in this case are: i) a = — 3r = — 6e, ii) a = — 2r = —2e, 
Hi) a = T. 

• Gauge invariant models. In this case we have r = and the only term remaining 
in Minkowski space-time is the gauge invariant one. The possibilities in this case 
are a = me with m = 0, —2, —4. 

Notice that except for the a = t = case (which is nothing but GR plus Maxwell 
electromagnetism) , all the considered cases break gauge invariance in general space- 
times. Therefore, there are six different classes of models which are indistinguishable 
from GR by means of Solar System experiments and, therefore, do not spoil the current 
bounds on the PPN parameters. To our knowledge best, none of these models (apart 
from Maxwell's theory) had been considered previously in the literature. 

We will also use the present constraints on the variation of the Newton's constant 
which are given by G/G ^ 10~^^ yr~^ [2]. 




(1) 



A^ = (AO,o,o) 



(2) 
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3. Classical and quantum stability 

To study the existence of unstable classical modes and ghosts we shall perform 
perturbations around a Minkowski background. In addition, we will also consider 
perturbations around the constant background vector field introduced in ([2]). This 
is possible because Minkowski space-time is an admissible solution of the theory in the 
case in which the vector field takes a constant value, as that introduced above. For this 
constant background, the vector field breaks Lorentz invariance and we have a preferred 
frame defined as that in which the vector field has only temporal component. In this 
background we decompose the perturbations in Fourier modes and solve the equations 
for the corresponding amplitudes. That way, we obtain the dispersion relation, which 
provides us with the propagation speed of the modes, that is required to be real in 
order not to have exponentially growing perturbations. As we explained above, here we 
shall not care about superluminal propagation of the modes, although it could be easily 
imposed at some point. Notice that although we are assuming constant background, in 
practice, this background could evolve on cosmological timescales. The effects of such 
evolution could have important cosmological consequences [6],[7j. 
On the other hand, we define the energy for the modes as [13]: 

"-(^"'-S^^S') (3) 

(2) (2) 

where T^jy and Gin, are the energy-momentum tensor of the vector field and the 
Einstein's tensor calculated up to quadratic terms in the perturbations and (■ ■ ■) 
denotes an average over spatial regions. Then, we insert the solutions obtained for the 
perturbations into this expression and study under which conditions they are positive. 

Notice also that the preferred frame respects the invariance under spatial rotations 
and therefore, in order to simplify the analysis we can perform the usual split of the 
perturbations into spin-0 (scalar), spin-1 (vector) and spin-2 (tensor). For simplicity, 
the longitudinal gauge has been chosen in the calculations below, although we have 
checked that the final results for the mode frequencies and energy densities defined in 
([3]) do not depend on the gauge choice. 

The scalar perturbations of the vector field can be written as 5"^ = (5*0, VS) and 
the perturbed metric in the longitudinal gauge is: 

ds'^ = (1 + 2(f))dt'^ - (1 - 2ij)6ijdx'dx^ (4) 

For vector perturbations we have = (0, v) and the metric is as follows: 

ds'^ = dt^ + 2F ■ dxdt - 6ijdx'dx^ (5) 

with V-'U = V--F = 0. Although the vector field does not generate tensor perturbations, 
we still can have effects by the presence of the background vector field as we shall see 
later. Let us first consider the scalar and vector perturbations. 
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3.1. Gauge non-invariant models 

3.1.1. Model I: a = — 3r = — 6e In this model, the Fourier components of the 
gravitational potentials relate to those of the vector field as follow^: 

AeA 



1 - 8eA^ L 
AeA 



3Sk — Si 



Ok 



[1 - 32eA^)Sk + S^ 



Ok 



1 - 8ev42 

whereas S^k and Sk are related to each other by means of: 



(6) 
(7) 



5, 



OA: 



A;2(l + 8eA2) 



ASeA^k'^ 



dSk 
dt 



[I - SeA^)(1 - 165^2) (1 - QAsA:' 



d^Sk 

) 

dt^ 



(8) 



Thus, the problem is solved once we know the solution of Sk, which happens to satisfy 
the following fourth order equation: 



d^S, 



+ 



dt^ ' (1- 166^2) (l-64ey42) 
This equation yields two independent modes: 

with their respective dispersion relations: 

oof = 



2(1 - AOeA' 



d^Sk 

) 

dt^ 



+ k^Sk 



0. 



(9) 
(10) 

(11) 
(12) 



2 l-lQeA^' 

Then, in order to have stable solutions we need to satisfy the following condition: 

UeA^ < 1 (13) 

because, in that case, both modes have real propagation speeds. 

On the other hand, the energy density evaluated over the solutions for each mode 
is given by: 



pis) 



l2^eA^) 



(1 - 64^^2)2 

l2 



\Ci\ 



(14) 
(15) 

(16) 

Now, by combining (JT3l) and (fT6l) we find the following viability condition for the scalar 
modes: 



These energies are both positive under the following constraint: 
1 < 128eA^ < 4 



1 < 128eA^ < 2 



(17) 



I Hereafter we will measure the field in units of vIttG. 
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Concerning vector perturbations, we obtain the following relation between the 
amplitudes: 



1 - 165^2'''= 

and both evolve as plane waves with the following dispersion relation: 



1 - 16eA'^ 
which leads to the stability condition: 

16eA^ < 1 

Finally, the energy density associated to the vector perturbations is given by: 



1- 326^2 2 

Se— T^r^k \vok\ 



[1 - IQey 

From this expression we see that s must be positive and we have to satisfy: 
32£ < 1 



(18) 

(19) 
(20) 
(21) 

(22) 



in order not to have vector modes with negative energy. 

In this model, according to the definition in the Appendix, G = 1 and the 
constraints on the variation of G do not set any limit on the possible variation of A. 



3.1.2. Model II: a — —2r — —2e In this case, and ijj are given in terms of the 
perturbation of the vector field as follows: 

3rA 



<t>k = 



1 - ?,tA 



"l-3rA2 

On the other hand, the perturbation 5*0^ can be expressed as: 



(23) 
(24) 



1 



2F(1 - 3rA2) 



^1-6tA')(1-15t^') 



d^Sk 



A;'(1 + 3t^') 



dSk 
dt 



(25) 



Then, all the perturbations are given in terms of Sk, for which we can obtain the 
following fourth-order differential equation: 



d^Sk 
dt^ 



+ 2k' 



1 - StA^ 



(l-6T>12)(i_i5rA2) 



2(1 - 9t^2 



' dt^ 



The solution of this equation is a superposition of two independent modes: 



—iU!-t 



0. (26) 
(27) 



with their respective dispersion relations: 



UJ. 



^ (1-6t^2)(i_i5^^2) 

±3|r|AV(l - 3tA2)(5 - 27rA2) 



[1 - 12tA^ + 27t2A^ 



(28) 
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Figure 1. This plot shows the dependence of both modes propagation speeds 
on tA^ for the model II: blue for lj^ and dashed-red for 



Then, we have two modes with two different speeds of propagation which depend on tA^ 
as it is shown in Fig. [H The a;_-mode has real propagation speed for tA^ < |, whereas 
for the u;+-mode to have real propagation speed we need to satisfy either tA'^ < ^ 
or tA"^ > |. Therefore, the necessary condition in order not to have instabilities is 
tA"^ < j^. Finally, the degeneracy disappears for tA"^ = when both propagation 
speeds are 1, recovering thus the usual result. 

The energy density corresponding to each mode can be expressed as: 

P^±^ = rfATA')k'\C^f (29) 

where /±(Ty4^) are the functions plotted in Fig. [2l Notice that / and r must have the 
same sign for the energy density to be positive. We find the following condition in order 
to have positive energy density for both modes: 

tA^ e (0.033, 0.105) U (0.383, 0.5) (30) 

For this model, the vector perturbation on the metric relates to that of the vector 
field by means of: 

which evolve as: 

Vk = vokc-"'^' (32) 
where vok ■ k = and 

From this expression we obtain that for tA"^ < | the propagation speed is real. 
The energy density corresponding to the vector perturbations is: 

P = (1 - 6rA^)2 l^o^l ^34) 



Viability of vector-tensor theories of gravity 



8 



f(TA-) 




Figure 2. This plot shows the functions /±(tA^) which determine the sign 
of the energy density for the scalar modes in the model II: blue for /-|_ and 
dashed-red for 

One can easily verify that this expression is positive if: 

tA^ G [0, c_] U [c+, oo) (35) 

with c± = |(1 ± ^). Note that r > is a necessary condition and that the singular 
value tA^ = I is not contained in the interval. 

D 

In this case, G = I + GtA"^. The present constraints on the variation of the Newton's 
constant will translate into a limit on the possible variation of A which will depend on 
the present cosmological value of A. 



3.1.3. Model III: a = t The perturbations of this model propagate at the speed of 
light so there are no classically unstable modes. However, the scalar perturbations are 
not just plane waves, but they have also a growing mode: 

Sok = [-ik\{Do + ikD)t + Do] e^*'* 

= [X{Do + ikD)t + D]e-'''' (36) 



with 



A 



e + 2r(2e + t)A'^ 



(37) 



e + r + 2r(2e + r)A2 

In principle, this solution can give rise to both positive or negative energies 
depending on the value of the amplitudes. However, there exists a way to get this 
difficulty around which is motivated by the fact that the action corresponding to this 
model can be rewritten as that of Maxwell's electromagnetism in the Gupta-Bleuler 
formalism, that is: 

1 



S 



-R - -F.^F"" + ^ (V.A"? 



(38) 



Notice that this action is the physically relevant action in the covariant formalism, 
since the energy of the modes and all the observables are calculated from (1381) (see for 
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instance [S]). In order to get rid of the negative energy modes, it is necessary to restrict 
the Hilbert space of the theory, by imposing the (Lorentz) condition {(j)\dfj,A^\(j)) = 0, 
which determines the physical states As shown in [7j, this is equivalent to remove 
the growing mode i.e. Dq = —ikD so that perturbations of the vector field propagate as 
pure plane waves. Then, using this condition we obtain that, as in the electromagnetic 
case, in the restricted Hilbert space, the energy of the scalar modes is identically zero. 
A detailed treatment of the quantization for this model is performed in [7j. 
The solution of the vector perturbation of the vector field is: 

Vk = vokc-'"' (39) 

Moreover, the vector perturbation of the metric vanishes. 

The energy density associated to the vector perturbations is: 

p(^) = (2£ + r)P|t;o,|' (40) 

which is positive if 25 + r > 0. 

In this model we have again G = 1 and no constraints on the variation of A can be 
established. 



3.2. Gauge invariant models: t = 0, a = me, m = 0, —2, —4 

As the case m = is nothing but Einstein's gravity plus electromagnetism (which 
can be easily seen to satisfy all the viability conditions), we shall focus just on the 
cases m = —2, —4. In such cases, we can obtain the following relations between the 
perturbations: 

■ _ 1 + (m + 3)meA^ 
^ ~ meA{l + 2meA^y^' 

_ l-meA^{3-2m^eA'^) 
°' ~ meA{l + 2meA^) 
(j)k =^k-2meASk. (41) 

Therefore, all the perturbations can be immediately obtained once we know the solution 
of ipk which happens to evolve as: 

ipk = Gue-'-'\ (42) 



with: 



= ii±^(i±I!!)!^,.. (43) 



' 3 1 + 2meA^ 
Then, the classical stability of these modes is guaranteed in the range: 



eA^ E -oo, U — — , oo . (44) 

V 2m J \ 2(4 + m)m' / ^ ' 

Notice that the second interval vanishes for the model with m = —4. 
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When we insert the corresponding solutions into the energy density we obtain: 

P ^ 27rG(l + 2meA-) ^ " ^^^^ 

This energy density is positive for: 

eA^ e (-00, a+) U (^-^, , (46) 

with: 

«± = 7 \ V (47) 

2m (l±^-^) 

Notice that a_ is 00 for the model m — —4. 

On the other hand, the vector perturbations of the field evolve as plane waves 
'^k — ■J^ofee""^"* with the following dispersion relation: 
„ l + (2 + ^m)m£A2 

which imposes the condition: 

sA' e (-OC, U 1 ^ , 00) (49) 

\ 2m / \ (2 + ^m)m J 

in order to have real propagation speed. Notice that for m = —4 the second interval 
vanishes. 

Besides, the vector perturbation of the metric relates to v by means of: 
- 2meA ^ 

;Vk. (50) 



1 + 2meA'^ 

These solutions for the vector perturbations lead to the following expression for the 
energy density: 

(.) 1 + ^msA^ (1 + (1 + \m)meA^) ^ . 

P ^ ' (l + 2m.A^) ^ l""'^" ^''^ 

The requirement for this energy to be positive is: 

1 
2m 
where 

^± = (0 ^ I ^ ^^3) 

m(2 ± \J —va) 

Note also that for m = — 4, fe_ becomes +00. 

In this case we get: G — \ — £m(4 ^rn)A?^ so the only model with G 1 is that 
with m = — 2. 



e (0,6+) u( -7^,6-1, (52) 
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Model I 


Model II 


Model III 


Gauge invariant models 
m = -2, -4 


Classical Stability 


GieA^ < 1 


15tA^ < 1 


Always 


sA^ ^ [ , 7 

>- 2m' (2+im)m 


Gravitational Waves 


2eA^ < 1 


12tA^ < 1 




eA^ < ^ 


Quantum Stability 


1 < 1286^2 < 4 


Not viable 


a^A^ = and 2e + r > 




Viability condition 


1 < 128eA'^ < 2 


Not viable 


a^A^ = and 2£ + r > 


Incompatible 



Table 1: In this table we summarize the conditions obtained in order to have both classical and 
quantum stability for the models with the same set of PPN parameters as GR studied in this work. 
The m = gauge invariant model satisfies all the viability conditions. 



4. Gravitational Waves 

At first glance, one may think that as the vector field does not generate tensor modes at 
first order, gravitational waves will not be affected. However, the presence of a constant 
value of the vector field in the background can modify the speed of propagation of tensor 
perturbations. For the general vector-tensor action ([T]) we have the following dispersion 
relation: 

= H ..2 (54) 

for both © and ® polarizations. Therefore, the speed of gravitational waves is modified 
in the presence of the vector field, recovering the usual value for A = 0. Thus, if 
2{(t — t)A^ > —1 we do not have unstable modes. In particular, the constraint a — T>0 
is a sufficient condition (although not necessary) which is independent of the background 
vector field. 

On the other hand, the energy density associated to the tensor perturbations is also 
modified by the presence of the background vector field: 

(,) _ l + 4(a-r)A^ 

P -32vrGl + 2(a-r)A2('^®' + ^ 

where C®, are the amplitudes of the corresponding graviton polarizations. 

Then, in order not to have modes with negative energy density we need either 
2((T — t)A'^ < —1 or 4(cr — t)A'^ > —1. These conditions combined with the classical 
stability condition lead to the constraint 4(cr — t)A'^ > —1. For models I and II this 
condition reads 2^74^ < 1 and VItA? < 1 respectively. On the other hand. Model III 
has 0" = r so gravitational waves are unaffected. Finally, for the gauge invariant models 
we obtain eA"^ < 
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5. Conclusions 

We have studied the viability conditions for unconstrained vector-tensor theories of 
gravity. In Table 1 we summarize the different conditions obtained for the different 
models analyzedl§| If we concentrate only on classical stability and local gravity 
constraints, we find that there are different types of vector-tensor theories which can be 
made viable for certain values of the coefficients. Thus for example: Model I with e < 
or Model II with r < for arbitrary values of A. 

However if we also impose that the models are free from ghosts, we have found 
that only theories of the Maxwell type or Maxwell plus a gauge-fixing term can be 
made compatible with all the consistency conditions for arbitrary A. Notice that on 
small (sub-Hubble) scales, such vector-tensor theories are indistinguishable from General 
Relativity, however it has been shown (see P, [7] for different examples) that in a 
cosmological scenario, the presence of the (dynamical) background vector field could 
have important consequences. Indeed in [7] it has been shown that the electromagnetic 
theory with a gauge fixing term can explain in a natural way the existence and the 
smallness of the cosmological constant, solving in this way the problem of establising 
what is the fundamental nature of dark energy. Present and forthcoming astrophysical 
and cosmological observations could thus help us determining what is the true theory 
for the gravitational interaction. 
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Appendix: PPN parameters for vector-tensor theories of gravity. 



The PPN parameters for the vector-tensor theory ([T]) are given by [2]: 
1 + AuA^ (1 + 2^^|±^) 



with: 
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/3 



e 

A 



1 - 4tuA2 (l - 

^(3 + 7) + ^e 



2e+T 
_8uj_\ 
2e+T) 

7(7-2)' 



1 + 



G 



ai = 4(1 - 7) [1 + (2e + r)A] + IGcuA'^Aa 



a2 = 3(1 - 7) 

"3 



l + -(2£ + r)A 



+ SioA^Aa - 2 



bA^ 
~G 



^ Ci = C2 = C3 = C4 = 

(l-4tjA2)(2e + a-2tj) 
(1 -4cuA2)(2£ + r) + 32cjM2 

1 

2A\a - r)2 - {2e + r) [1 - AA'^^uj + a ^ r)] 
{2e + r)(l - 37) + 2(a - r)(l - 27) 



(2cu + a - r) [(27 - 1)(7 + 1) + 6(7 - 2)] 
-(27 -1)2(2.; + a) (1-^) r^O 



r = 

Moreover, it is possible to define an effective Newton's constant given by: 

^eff = G 

- 2A\a-r)il + e)]~\ 
The above expressions are obtained assuming Gg// = 1- 



-(7 + 1) + 6cuA2(7-1) 



